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^3 ■ ABSTRACT 
<. 

The state of the art of the theoretical calculations for heavy quarks photoproduc- 
' tion is reviewed. The full fixed order next-to-leading order massive calculation 

and the resummation of large \og(j)T /m) terms for differential cross sections are 
^ ' described. The implementation of a non-perturbative fragmentation function 

, describing the c—^D* meson transition is also discussed. 

00 

(yQ • 1. Introduction 

o ■ 



Heavy quarks production processes provide a powerful insight into our under- 
standing of Quantum Chromodinamics. The large mass of the heavy quark can make 
^ the perturbative calculations reliable, even for total cross sections, by cutting off in- 

p • frared singularities and by setting a large scale at which the strong coupling can be 

^ ■ evaluated and found - possibly - small enough. On the experimental side, the possi- 

bility to tag heavy flavoured hadrons by means of microvertex detectors can on the 
other hand provide accurate measurements. 
^ \ All these potentialities must of course be matched by accurate enough theoretical 

\ evaluations of the production cross section. In this talk I shall describe the state of 

the art of such calculations for heavy quarks photoproduction. I shall first review 
the next-to-leading order (NLO) QCD evaluations recently presented by Frixione, 
Mangano, Nason and Ridolfi. These calculations, available for total cross sections, 
one-particle and two-particles distributions, are now a consolidated result and provide 
a benchmark for future developments. 

Large logarithms appear in the NLO fixed order calculations and potentially make 
it less reliable in some regimes: log(S'/m^) and log(p^/m^) become large when the 
center of mass energy or the transverse momentum of the observed quark is 
much larger than its mass. I shall describe the resummation of log(p|./m^) terms, 
leaving the high energy resummation to Marcello Ciafaloni's talkB. 

The perturbative fragmentation function technique used in the resummation of 
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the large log(p|^/m^) terms has a non-perturbative extension which can be used to 
describe the transition from c quarks to D* mesons. I shall therefore also discuss the 
determination of these non-perturbative fragmentation functions and their inclusion 
into the heavy quarks photoproduction calculation, showing a comparison with data 
from HERA. 

2. Fixed Order NLO Calculation 

Heavy quarks photoproduction at leading order in the strong coupling looks a 
very simple process: only the tree level diagram 7(7 — QQ contributes at the partonic 
level, and the final answer for the total cross section is simple and well behaved, being 
finite everywhere. 

At a deeper thinking, however, problems seem to arise. For instance, one may ask 
himself why not to include initial state heavy quarks, coming from the hadron and to 
be scattered by the photon, like 7Q — > Qg. To include consistently such a diagram 
is not an easy task, especially if one wants to keep the quark massive. Taking it 
massless, on the other hand, would not only be a bad approximation but would also 
produce a divergent total cross section. 

A way out of this problem was provided by Collins, Soper and Stermani, who ar- 
gued that the following factorization formula holds for heavy quarks hadropro duct ion 
total cross sections: 

(y{^, m)=J2j fi/Hjj/H2^iij QQ; v^, m). (1) 

The sum on the partons runs only on i and j being gluons or light quarks, and 
the heavy quarks are only generated at the perturbative level by gluon splitting. 
There is therefore no need to try to accommodate them in the colliding hadrons 




and the relevant kinematics can be kept exact. Eq. (|I]) provides the basis for an 
exact perturbative calculation of heavy quarks production to NLO. For what concerns 
photoproduction, such a calculation has been first performed by P. Nason and K. Ellis, 
and subsequently confirmed by J. Smith and W.L. van Neerven i. 

When going to order aa^ in photon-hadron collision, however, a new feature 
appears. The photon can now couple directly to massless quarks, for instance in 
processes like 7g QQq, and in a given region of phase space a coUinear singularity 
will appear. It can be consistently factored out, but this requires the introduction 
of photon parton distribution functions (PDF) which, pretty much like the hadron 
ones, will describe the probability that before the interaction the photon splits into 
hadronic components (light quarks or gluons, in this case). Such a behaviour is 
sometimes called resolved photon (as opposed to direct). A full NLO calculation 
for heavy quarkphotoproduction will therefore also require a NLO calculation for 
hadroproductionE, where one of the PDF's will be the photon's one. A factorization 
scale n^, related to the subtraction of the singularity at the photon vertex, will link 
the two pieces and its dependence on the result will only cancel when both are taken 
into account. 

Frixione, Mangano, Nason and Ridolfii (FMNR) have recently presented Mon- 
tecarlo integrators for these two calculations, thereby allowing detailed comparisons 
with experimental data. A very extensive collection of such comparisons is presented 
in a recent revie J, from which we select some plots to be shown here. 

A comparison of total cross section experimental results and theoretical predic- 
tions for cc photoproduction is shown in fig. |l|. Although large uncertainties are 
present, the comparison suggests agreement between theory and experiment. The 
new HERA data, at large center of mass energy, can be seen to appear larger than 
the pointlike (= direct) photon prediction only. This suggests the need for a resolved 
photon component, but by no means can determine it precisely. 

One-particle transverse momentum (pt) distributions are shown in fig. The 
pure QCD predictions can be seen to be significantly harder than the data. However, 
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Fig. 3. Two particles correlations in fixed target experiment^. 



when corrected with two non-perturbative contributions they can be matched to 
the data. These non-perturbative addictions are meant to represent a primordial 
transverse momentum fc^ of the colhding partons, other than the one aheady taken 
into account by the QCD radiative corrections, and the effect of the fragmentation of 
the produced heavy quark into the observed heavy flavoured hadrons, here described 
by the so-called Peterson fragmentation function with e = 0.06. 

Comparisons between data and theory for two-particle correlations, like the az- 
imuthal difference A0 or the relative transverse momentum Pt{QQ) of the produced 
heavy quark pair, are shown in fig. ^ Distributions like these are trivial in lead- 
ing order QCD, since the Q and the Q are produced back-to-back. Hence, A0 = n 
and Pt{QQ) = 0. NLO corrections (as well as non-perturbative contributions) can 
broaden these distributions, and one could think of being able to perform a direct 
measurement of 0{al) effects. The plots do however show that non perturbative 
contributions play a key role in allowing a good description of the data. One can, 
however, still check that the same inputs allow for a good description of both one- 
and two-particles distributions, as seems to be the case here. 

The overall result of these comparisons can therefore be summarized as follows. 
Total cross sections seem to be well reproduced by the calculation both at fixed target 
and HERA regimes, but the huge uncertainties present both on the experimental and 
the theoretical side do not allow the study of finer details like, for instance, the 
determination of the resolved component at HERA. For what concerns transverse 
momentum distributions at fixed target, they can be reproduced after allowing for 
heavy quark fragmentation effects and for a primordial transverse momentum of the 
incoming partons of the order of 1 GeV. These same non-perturbative corrections 
also allow for a description of two-particles correlations, thereby pointing towards a 
consistent picture. 
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Fig. 4. Comparison between_fixed order (FMNR) and resummed (PFF) calculation for charm pho- 
toproduction pT distributioru. It is worth noticing how the two calculations describe differently the 
(unphysical) resolved and direct components, but agree on their sum (a physical observable). 



3. Large Transverse Momentum Resummation 

Like any perturbative expansion, the NLO calculation for heavy quarks photo- 
production is only reliable and accurate as long as the coefficients of the coupling 
constant remain small. Large terms of the kind log(p^/m^) do however appear in the 
cross section, and for growing pj- they will eventually became large enough to spoil 
the convergence of the series. Such terms need therefore to be resummed to all orders 
to allow for a sensible phenomenological prediction. Such a resummation has been 
performed along the following linesi. 

One observes that in the large-pr hmit {pt ^ ^) the only important mass terms 
are those appearing in the logs, all the others being power suppressed. This means 
that an alternative description of heavy quark production can be achieved by us- 
ing massless quarks and providing at the same time perturbative distribution and 



fragmentation functions also for the heavy quark, describing the logarithmic mass 
dependence. The factorization formula becomes 

da{pT) = YlJ Fi/HiifJ', [m])Fj/H2(/i, [Tn])da{ij k;pT, IJ,)D'^ {fi,m), (2) 

ijk 

with parton indices i,j and k also running on Q, taken massless in a, now an MS 
subtracted cross section for light partons production. The dependence on m of the 
parton distribution functions Fi/n, shown among square brackets in eq. (0), is only 
there when i or j happens to be the heavy quark Q. 

The key point is that the large mass of the heavy quark allows for the evaluation 
in perturbative QCD (pQCD) of its distribution and fragmentation functions. Initial 
state conditions for Fq/ni^o = and -D^ (/xo — can be calculated in pQCD at 
NLO level in the MS scheme: 

Fqih{x,^io = m) = (3) 
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(4) 



^?(^, /^o) = + (1 - x)^) log 4 (5) 

^S.-o(^""o) = (6) 

The massive logs will hence appear only through these function, which can then be 
evolved with the Altarelli-Parisi equations up to the large scale set by yU ~ pt- This 
evolution will resum to all orders the large logarithms previously mentioned. 

It is important to mention that due to the neglecting of power suppressed mass 
terms this approach becomes unreliable when ~ m. In this region only a case 
by case comparison with the full NLO massive calculation - here reliable and to be 
taken as a benchmark - can tell how accurate the resummed result is. 

Phenomenological analyses show that the effect of the resummation becomes size- 
able only at very large pt, say greater than 20 GeV for charm photoproduction. 
Fig. ^ shows the effect of such a resummation for a fixed photon energy in HERA-like 
kinematics. The resummed calculation can be seen to match the fixed order one at 

~ m, where resummation effects are not expected to be important, and to be- 
have more softly in the large region. This particular theoretical refinement should 
therefore not be phenomenologically overly relevant for present-day HERA physics, 
data being only available up to pt — 12 GeV. 



4. On the Inclusion oi c ^ D* Fragmentation Effects 

When comparing theory with data, one always faces the problem of describing 
as closely as possible what the experiments do observe. With heavy quarks pro- 
duction the problem lies in the experiments actually seeing the decay products of 



heavy flavoured hadrons rather than the heavy quark itself. This is due to the heavy 
quark strong and non-perturbative binding into a hadron prior to decay. This bind- 
ing involves the exchange and radiation of low-momentum (order Kqcd) gluons, and 
typically degrades the momentum of the hadron with respect to the one of the origi- 
nal quark. Such a degradation can be described with the help of a non-perturbative 
fragmentation function (FF) which, lacking the theoretical tools to calculate, can be 
extract by fitting experimental data. 

An often employed parametrization for such a function is the so called Peterson0 
one, which reads 

Dnp{z]e) — -2. (7) 

z[l-l/z-e/{l- z)] 

The value of e isDredicted to scale like Ag^^^/m^. For charm to D* fragmentation 
a global analysis" based on leading order Montecarlo simulations gives the value 
e ^ 0.06. This value has so far usually been taken as the reference one, and used for 
instance together with the NLO fixed order calculation by FMNR in the plots shown 
in Section 2. 

One should however carefully consider how e has been extracted from e^e^ ex- 
perimental data. Experiments usually report the energy or momentum fraction [xe 
or Xp) of the observed hadron with respect to the beam energy. On the other hand 
the fraction which appears as the argument of the non-perturbative FF is rather to 
be taken with respect to the fragmenting quark momentum, usually denoted by z 
(see for instance til for a discussion on this point). These two fractions are not co- 
incident, due to hard radiation processes which lower the momentum of the quark 
before it fragments into the hadron. In order to deconvolute these effects one usually 
runs a Montecarlo simulation of the collision process at hand, including both the per- 
turbative parton showers and the subsequent hadronization of the partons into the 
observed hadrons. The latter can be parametrized in the Montecarlo by the Peterson 
fragmentation function, and the value of e which best describes the data can be ex- 
tracted. Clearly this procedure leads to a resulting value for e which depends on the 
details of the description of the perturbative part. Indeed, the showering softens the 
momentum distribution of the heavy quark, producing an effect qualitatively similar 
to that of the non-perturbative FF. On the quantitative level, the amount of softening 
(and hence the value of e) required by the non-perturbative FF to describe the data 
is related to the amount of softening already performed at the perturbative level. A 
leading or a next-to-leading description of the showering can therefore produce differ- 
ent values for e, whose value is then not a "unique" and "true" one, but rather closely 
interconnected with the details of the description of the pQCD part of the problem. 

In ref. @ fits to D* data taken by the ARGUS and OPAL experiments have been 
performed with NLO accuracy using a fragmentation description for the heavy quark 
production like the one described in Section 3, complemented with the inclusion of a 




0. 1 0.2 IJ.3 0.4 (1.5 (1.6 (».7 O.S 0.9 1 0. 1 0.2 0.3 (1.4 0.5 0.6 0.7 0.8 0.9 



Fig. 5. Distributions of Djl-imesons as measured by the ARGUS and OPAL experiments, together 
with the theoretical curvest3 fitted to the same data with the (1 — x)°'x^ (full line) and the Peterson 
(dashed line) non-perturbative fragmentation functions. 

non-perturbative component via the ansatz 

D^\^,) = Dli^^) ® (8) 

represented by the convolution of a perturbatively calculable fragmentation function 
of the parton k into the heavy quark c and the non-perturbative form Dj^ describing 
the c —>■ D* transition. This non perturbative form is taken to be scale independent, 
i.e. all scaling effects are assumed to be described by the Altarelli-Parisi evolution of 
the perturbative part DKfi). A similar approach had already been introduced in §. 

Results for these fits are shown in fig. ^. The value for e has been consistently 
found to be of order 0.02 rather than the customary 0.06 one, resulting instead from 
fits with leading order evolution. Recalling the previous discussion, this comes to 
no surprise: next-to-leading order evolution softens more the heavy quark spectrum, 
and a harder non-perturbative fragmentation function is therefore needed to provide 
a satisfactory description of the data (see El for a full discussion). 

Similar fits to e~^e~ data have also been performed by Binnewies, Kniehl and 
KramerEl (BKK). These authors do instead find, again with NLO evolution, a value 
for e still close to the usual 0.06. This discrepancy, beyond irrelevant nomenclature 
differences, can be traced back to a discrepancy in the implementation of the factor- 
ization scheme. The scheme used in 0, as originally set up in i, is the customary 
MS one. Considering for instance the dominant non-singlet component only for sim- 
plicity, the e~^e~ QX momentum distribution da/dx can be schematically written 
as the convolution (= product in Mellin moments space) of a short distance coef- 
ficient function, an Altarelli-Parisi evolution kernel i?(/i, /xq), a perturbative initial 
state condition for the heavy quark perturbative fragmentation function (PFF) and 
a fixed non-perturbative FF, 

da{y/S, m) = (l + as{i^)c{\^ , j2)jE{fi, jjo)(j- + as{l^o)d{l^o,rn)^Dnp, (9) 



where the perturbative expansions of the coefficient function and the PFF have been 
exphcitly shown. The factorization scale /i is taken of the order of the (large) collision 
energy "s/S*, and the initial scale /xq is taken of the order of the quark mass m. 

BKK on the other hand, employing a scheme introduced by Kniehl, Kramer and 
Spira0 (KKS), write da{y/S,m) as 

da{y/S, m) = (l + as{ij)c{y/S, fi) + a^(/i)(i(/io, m)^E{^, f^o)Dnp- (10) 

These two expressions can be seen to differ by O(a^) terms. However, one of these 
terms is given by 

c^sM - as(/io) = -&o«s log — (11) 

and is, therefore, one of the next-to-leading logarithms (NLL) a'^ \og''~^ {VS / m) we 
are resumming. Hence the two calculations differ by a NLL term and cannot possibly 
both implement correctly a resummation at the NLL level. 

To better understand the discrepancy, the BKKS scheme can for instance be 
rewritten in the form (P), with an initial state condition for the PFF containing the 
large scale fi as the argument for as rather than the small one /xq- This choice of 
a large scale is however in contradiction with the factorization theorem hypotheses, 
which only allow for small scales in initial conditions, to avoid the appearance of 
unresummed large logs. Choosing the large /i leads at a practical level to the differ- 
ence being reabsorbed into a different value for the e parameter, which happens quite 
accidentally to be 0.06 rather than 0.02. One can show that, replacing in the BKKS 
formula (|ToD the as{fi)d{fio,m) term with as{iiQ)d{no,m) (or alternatively appropri- 
ately modifying the NLO splitting vertices in the evolution kernel), e = 0.02 is once 
again found from the NLL fits within this scheme too. 

On the phenomenological side, and making use of the universality argument, one 
can now argue that the use of a "harder" Peterson form with e = 0.02 is probably 
more suited when combined with a NLO perturbative calculation like the FMNR 
one which, albeit only at fixed order, contains NLL gluon radiation. Decreasing 
e means increasing the cross section at large pt, being the pt distribution steeply 
falling with increasing transverse momentum. This could help reconciling the HERA 
experimental dataS with the perturbative NLO calculation, which was shown to 
underestimate them a little when convoluted with a Peterson with e = 0.06: fig. ^ 
shows, on the left, how the cross section for D* photoproduction at HERA increases 
with decreasing e and, on the right, a comparison of the HI data with the fixed order 
prediction by FMNR (e = 0.06) and the fragmentation functions one with e = 0.02. 
One should notice that the pt values involved are still pretty small: this means that 
the fixed order calculation is still reliable and the accuracy of the resummed one has 
to be assessed ffist by comparing with the former. In this case they are found to be 
in good agreement, the difference in the plot being mainly given by the different e 
values. 
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Fig. 6. Effect on the D* photoproduction cross section of a decreasing value for e (left), and 
comparison of HI data with the fixed order calculation by FMNR (histograms, e — 0.06) and the 
fragmentation function approach (smooth lines, e = 0.02). 



Last but not least, it is worth mentioning how, going from LO to NLO analy- 
ses, a similar hardening of the non-perturbative fragmentation function is also ex- 
pected for the b quark. The corresponding increase of the hadroproduction bottom 
Pt distributionsS would be welcome in the light of the Tevatron data presently over- 
shooting the theoretical predictions by at least 30%. 
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